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Abstract 

Two planar graphs G\ and Gi sharing some vertices and edges are simultaneously planar if they 
have planar drawings such that a shared vertex [edge] is represented by the same point [curve] in both 
drawings. It is an open problem whether simultaneous planarity can be tested efficiently. We give a 
linear-time algorithm to test simultaneous planarity when the two graphs share a 2-connected subgraph. 
Our algorithm extends to the case of k planar graphs where each vertex [edge] is either common to all 
graphs or belongs to exactly one of them, and the common subgraph is 2-connected. 
Keywords: Simultaneous Embedding, Planar Graph, PQ Tree, Graph Drawing 



1 Introduction 

Let G\ = (Vi,Ei) and G2 = (^,£2) be two graphs sharing some vertices and edges. The simultaneous 
planar embedding problem asks whether there exist planar embeddings for G\ and G2 such that, in the 
two embeddings, each vertex v G V\ n V2 is mapped to the same point and each edge e G E\ n E2 is 
mapped to the same curve. We show that this problem can be solved efficiently when the common graph 
(Vi n V 2 , Ei n E 2 ) is 2-connected. 



The study of planar graphs has a long history and has generated many deep results 30-32 . There 
is hope that some of the structure of planarity may carry over to simultaneous planarity. A possible 
analogy is with matroids, where optimization results carry over from one matroid to the intersection of two 
matroids [12]. On a more practical note, simultaneous planar embeddings are valuable for visualization 
purposes when two related graphs need to be displayed. For example, the two graphs may represent 
different relationships on the same node set, or they may be the "before" and "after" versions of a graph 
that has changed over time. 

Over the last few years there has been a lot of work on simultaneous planar embeddings [2|[3| |13[[l4| 



17-20 27J, often under the name "simultaneous embeddings with fixed edges". We mention a few results 
here and give a more detailed description in the Background section below. A major open question is 
whether simultaneous planarity of two graphs can be tested in polynomial time. The problem seems to be 
right on the feasibility boundary. The problem is NP-complete for three graphs (20 and the version where 
the planar drawings are required to be straight-line is already NP-hard for two graphs and only known 
to lie in PSPACE |15|. On the other hand, several classes of (pairs of) graphs are known to always have 



simultaneous planar embeddings 113 14 18 ,19, 271 and there are efficient algorithms to test simultaneous 



planarity for some very restricted graph-classes: biconnected outerplanar graphs 1181 , the case where one 
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graph has at most one cycle JTT] or the case where the common graph is a star |3j. The last result me 
mention here is that the problem of whether a planar embedding of a subgraph can be extended to a planar 
embedding of the whole graph is equivalent to testing simultaneous planarity between a planar graph and 
a 3-connected graph, and a linear time algorithm for these problems was given by Angelini et al. 12]. 

In this paper we give a linear time algorithm to test simultaneous planarity of any two graphs that 
share a 2-connected subgraph. Our algorithm builds on the planarity testing algorithm of Haeupler and 
Tar j an |22| , which in turn unifies the planarity testing algorithms of Lempel-Even-Cederbaum [29], Shih- 
Hsu |34| and Boyer-Myrvold [9] . We show how to extend our algorithm to the case of k graphs where each 
vertex [edge] is either common to all graphs or belongs to exactly one of them, and the common graph 
is 2-connected. Although we concentrate on the decision version of the problem, we also show that two 
simultaneous planar graphs have a simultaneous planar drawing in which one graph is straight-line planar, 
and each edge of the other graph is drawn as a polygonal line with at most \Vx D V2I bends. 

Independently (and at the same time as our conference version (2l]), Angelini et al. [3] gave an efficient 
algorithm to test simultaneous planarity of any two graphs that share a 2-connected subgraph. Their 
algorithm is based on SPQR-trees, takes linear time (in the final version [4]), and is restricted to the case 
where the two graphs have the same vertex set. Very recently Blasius and Rutter [6] have generalized our 
idea of simultaneous PQ trees, and shown how to test simultaneous planarity of two biconnected graphs 
that share a connected subgraph. This is a substantial strengthening of our result. 

The paper is organized as follows: Section [l.l| gives more background and related work. In Section [2] 
we review and develop some techniques for PQ-trees, which are needed for our algorithm, and in Section [3] 
we review the Haeupler- Tarj an planarity testing algorithm. Section [4] contains our simultaneous planarity 
testing algorithm, including the extension to k graphs, and the result on drawing two simultaneous planar 
graphs. 



1.1 Background 

Versions of simultaneous planarity have received much attention in recent years. Brass et al. [10] intro- 
duced the concept of simultaneous geometric embeddings of a pair of graphs — these are planar straight-line 
drawings such that any common vertex is represented by the same point. Note that a common edge will 
necessarily be represented by the same line segment. It is NP-hard to test if two graphs have simulta- 
neous geometric embeddings 15 1. For other work on simultaneous geometric embeddings see [I] and its 



references. 

The generalization to planar drawings where edges are not necessarily drawn as straight line segments, 
but any common edge must be represented by the same curve was introduced by Erten and Kobourov |14| 
and called simultaneous embedding with consistent edges. Most other papers follow the conference version 
of Erten and Kobourov's paper and use the term simultaneous embedding with fixed edges (SEFE). In our 
paper we use the more self-explanatory term "simultaneous planar embeddings." A further justification 
for this nomenclature is that there are combinatorial conditions on a pair of planar embeddings that are 
equivalent to simultaneous planarity. Specifically, Jiinger and Schultz give a characterization in terms 
of "compatible embeddings" [Theorem 4 in [27]]. Specialized to the case where the common graph is 
connected, their result says that two planar embeddings are simultaneously planar if and only if the cyclic 
orderings of common edges around common vertices are the same in both embeddings. 



Several papers 13 14 19 show that pairs of graphs from certain restricted classes always have simulta- 
neous planar embeddings, the most general result being that any planar graph has a simultaneous planar 
embedding with any tree Il9j. On the other hand, there is an example of two outerplanar graphs that 
have no simultaneous planar embedding |19|. The graphs that have simultaneous planar embeddings when 
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paired with any other planar graph have been characterized 18 . In addition, Jiinger and Schultz |27 
characterize the common graphs that permit simultaneous planar embeddings no matter what pairs of 
planar graphs they occur in. 

There are efficient algorithms to test simultaneous planarity for pairs of biconnected outerplanar 
graphs [18], and for pairs consisting of a planar graph and a graph that is either 3-connected [2] or 
has at most one cycle [IT] . Testing simultaneous planarity when the common graph is a tree was shown 
to be equivalent to a book-embedding problem which together with |23| implies that testing simultaneous 
planarity when the common graph is a star can be done in linear time |3j. 

There is another, even weaker form of simultaneous planarity, where common vertices must be rep- 
resented by common points, but the planar drawings are otherwise completely independent, with edges 
drawn as Jordan curves. Any set of planar graphs can be represented this way by virtue of the result that 
a planar graph can be drawn with any fixed vertex locations |33| . 

The idea of "simultaneous graph representations" has also been applied to intersection representations. 
For the case of interval graphs, the problem is to find representations for two graphs as intersection graphs 
of intervals in the real line, with the property that a common vertex is represented by the same interval in 
both representations. This can be done in polynomial time |25|. There are also polynomial time algorithms 



to find simultaneous representations of pairs of chordal, comparability and permutation graphs |26|. 



2 PQ-trees 

Many planarity testing algorithms in the literature |9, 22, 29, 34] use PQ-trees (or a variation) to obtain 
a linear-time implementation. PQ-trees were discovered by Booth and Lueker [8] and are used, not only 
for planarity testing, but for many other applications like recognizing interval graphs or testing matrices 



for the consecutive-ones property. We first review PQ-trees and then in Subsection 2.1 we show how to 
manipulate pairs of PQ-trees. 

A PQ-tree represents the permutations of a set of elements satisfying a family of constraints. Each 
constraint specifies that a certain subset of elements must appear consecutively in any permutation. PQ- 
trees are rooted trees with internal nodes being labeled 'P' or 'Q', and are drawn using a circle or a 
rectangle (or double circles in [24]), respectively. The leaves of a PQ-tree correspond to the elements 
whose orders are represented. PQ-trees are equivalent under arbitrary reorderings of the children of a 
P-node and reversals of the order of children of a Q-node. We consider a node with two children to be 
a Q-node. A leaf-order of a PQ-tree is the order in which its leaves are visited in an in-order traversal 
of the tree. The set of permutations represented by a PQ-tree is the set of leaf-orders of equivalent PQ- 
trees. Given a PQ-tree tree T on a set U of elements, adding a consecutivity constraint on a set S C U, 
reduces T to a PQ-tree T", such that the leaf-orders of T' are precisely the leaf-orders of T in which the 
elements of S appear consecutively. Booth and Lueker [8] gave an efficient implementation of PQ-trees that 
supports this operation in amortized Od^l) time. Their implementation furthermore allows an efficient 
complement-reduction, i.e., adding the constraint on a set S when U \ S is given in amortized 0(\U \ S\) 
time 
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Although PQ-trees were invented to represent linear orders, they can be reinterpreted to represent 



circular orders as well 22 : Given a PQ-tree we imagine that there is a new special leaf s attached as the 
"parent" of the root. A circular leaf order of the augmented tree is a circular order that begins at the 
special leaf, followed by a linear order of the remaining PQ-tree and ending at the special leaf. Again a 
PQ-tree represents all circular leaf-orders of equivalent PQ-trees. It is easy to see that a consecutivity 
constraint on such a set of circular orders directly corresponds to a consecutivity constraint on the original 
set of linear leaf-orders and a reduction on the circular orders corresponds to a standard or complement 
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reduction. Note that using PQ-trees for circular orders requires solely this different view on PQ-trees but 
does not need any change in their implementation. As such it turns out that PC-trees introduced in |24| 
are the exact same data structure as PQ-trees albeit with this circular interpretation. 

2.1 Intersection and projection of PQ-trees 

In this section we develop simple techniques to obtain consistent orders from two PQ-trees. More precisely 
when two PQ-trees share some but not (necessarily) all leaves, we want to find a permutation represented 
by each of them with a consistent ordering on the shared leaves. 

Note that the set of such consistent orderings is not representable by a PQ-tree. For example, consider 
the ground set {1,2,3,4} and the constraint that {2,3} must be consecutive among {1,2,3}. This is 
why we restrict our goal to finding one permutation of each PQ-tree such that the two permutations are 
consistent on the shared leaves. 

The idea is to first project both PQ-trees to the common elements, intersect the resulting PQ-trees, 
pick one remaining order and finally "lift" this order back. We now describe the individual steps of this 
process in more detail. 

The projection of a PQ-tree on a subset of its leaves S is a PQ-tree obtained by deleting all elements not 
in S and simplifying the resulting tree. Simplifying a tree means that we (recursively) delete any internal 
node that has no children, and delete any node that has a single child by making the child's grandparent 
become its parent. This can easily be implemented in linear time. 

Given two PQ-trees on the same set of leaves (elements) we define their intersection to be the PQ-tree 
T that represents exactly all orders that are leaf-orders in both trees. This intersection can be computed 
in polynomial time as follows. 

1. Initialize T to be the first PQ-tree. 

2. For each P-node in the second PQ-tree, reduce T by adding a consecutivity constraint on all the 
P-node's descendant leaves. 

3. For each Q-node in the second tree, and for each pair of adjacent children of it, reduce T by adding 
a consecutivity constraint on all the descendant leaves of the two children. 

Using the efficient PQ-tree implementation, computing such an intersection be sped up to amortized 
linear time in the size of the two PQ-trees. This is a relatively straight-forward extension of the PQ-tree 
reduction and described in detail in Booth's thesis [7]. 

These two operations are enough to achieve our goal. Given two PQ-trees T\ and T2 defined on different 
element (leaf) sets, we define S to be the set of common elements. Now we first construct the projections 
of both PQ-trees on S and then compute their intersection T as described above. Any permutation of S 
represented by T can now easily be "lifted" back to permutations of T\ and of T2 that respect the chosen 
ordering of S. Furthermore, any two permutations of T\ and T2 that are consistent on S can be obtained 
this way. 

We note that techniques to "merge" PQ-trees were also presented by Jiinger and Leipert [28] in work 
on level planarity. Their merge is conceptually and technically different from ours in that the result of their 
merge is a single PQ-tree whereas we compute two orderings that are consistent on common elements. The 
set of all these consistent orderings cannot be captured by a PQ-tree. 
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3 Planarity testing 



In this Section, we review the recent algorithm of Haeupler and Tarjan [22] for testing the planarity of a 
graph. Next we extend it to an algorithm for testing simultaneous planarity. We first begin with some 
basic definitions. 

Let G = (V, E) be a connected graph on vertex set V = {v%, ■ ■ ■ , v n } and let O be an ordering of the 
vertices of V. An edge ViVj is an in-edge of Vi (in O) if Vj appears before Vi in O, and ViVj is an out-edge 
of Vi if Vj appears after Vi in O. 

An st-ordering of G is an ordering O of the vertices of G, such that the first vertex of O is adjacent to 
the last vertex of O and every intermediate vertex has an in-edge and an out-edge. It is well-known that 
G has an st-ordering if and only if it is 2-connected. Further, an st-ordering can be computed in linear 
time 
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The planarity testing algorithm of Haeupler and Tarjan embeds vertices (and their adjacent edges) one 
at a time and maintains all possible partial embeddings of the embedded subgraph at each stage. Edges 
with two, one or no endpoint(s) embedded are called (fully-)embedded, half-embedded and non-embedded 
respectively. For the correctness of the algorithm it is crucial that the vertices are added in a leaf-to-root 
order of a spanning-tree. This guarantees that, at any time, the remaining vertices induce a connected 
graph and hence lie in a single face of the partial embedding which can be thought of as the outer face. 
The crucial observation is that now (for the sake of extending this embedding) a partial embedding can 
be completely characterized by describing, for each connected component, the cyclic order of the half- 
embedded edges around it. Furthermore, PQ-trees can be used in a natural way to describe the set of all 
possible cyclic orders and thus all possible partial planar embeddings. Before describing how to update the 
PQ-trees corresponding to sets of partial embeddings when a vertex is added, we first discuss the choice of 
the specific embedding order we use. In general using either a leaf-to-root order of a depth-first spanning 
tree or an st-order leads to particularly simple implementations that run in linear-time. Indeed these two 
orders are essentially the only two orders in which the algorithm runs in linear-time using the standard 
PQ-tree implementation. Our simultaneous planarity algorithm will use a mixture of the two orders: We 
first add the vertices that are contained in only one of the graphs using a depth-first search order and 
then add the common vertices using an st-ordering. Note that if all vertices are common, both G\ and G2 
will be 2-connected and have a common st-order. In this case we only use this st-order. We now give an 
overview of how the update steps of the planarity test work for each of these orderings. 



Leaf-to-root order of a depth-first spanning tree: 

Let vi, i>2, • • ■ , v n be a leaf-to-root order of a depth-first spanning tree of G. Note that at stage i, the vertices 
{v 1, • • • , Vi} may induce several components. We maintain a PQ-tree for each component representing the 
set of possible circular orderings of its out-edges. Using a depth-first spanning tree, in contrast to an 
arbitrary spanning tree, has the advantage that we can easily maintain the invariant that the edge to the 
smallest node greater than i will be the special leaf. 

Adding Vi+i can lead to merging several components into one. It is easy to see, e.g., in the left part 
of Figure [2j that after embedding Uj+i the edges of each merging component that go to Vi + \ have to be 
consecutive since otherwise a half-embedded edge would be enclosed in the newly formed component. To 
go to the next stage, we thus first reduce each PQ-tree corresponding to such a component by adding a 
consecutivity constraint that requires the set of out-edges that are incident to Vi + \ to be consecutive. We 
then delete these edges since they are now fully embedded. By the invariant stated above the special leaf is 
among these edges. Note that the resulting PQ-tree for a component now represents the set of all possible 
linear-orders of the out-edges that are not incident to Ui+i. Now we construct the PQ-tree for the new 
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merged component including Vi+i as follows: 

Let vi be the parent of Vi+\ in the depth-first spanning tree. The PQ-tree for the new component 
consists of the edge Vi + \Vi as the special leaf and a new P-node as a root and whose children are all 
the remaining out-edges of Vi+i and the roots of the PQ-trees of the reduced components (similar to the 
picture in Figure [I]). Note that by choosing the edge Ui+iUj as the special leaf we again maintain the above 
mentioned invariant. 

It is easy to verify that these operations capture exactly all possible circular orders of half-embedded 
edges around the new embedded connected component. Note that if the reduction step fails at any stage 
then the graph must be non-planar. Otherwise the algorithm concludes that the graph is planar. 

st-order: 

Let vi,V2,--- ,v n be an st-order of G. This order is characterized by the fact that at any stage % 6 
{1, • • • ,n — 1} not just the non-embedded vertices but also the embedded vertices {v±,--- ,Vi} induce a 
connected component. This results in the algorithm having to maintain only one PQ-tree T, to capture all 
possible circular orderings of out-edges around this one component. Furthermore since v\v n is an out-edge 
at every stage, it can stay as the special leaf of Tj for all %. 

The update step for embedding the next vertex is the same as above. At stage 1, the tree T\ consists 
of the special leaf v\v n and a P-node whose children are all other out-edges of v±. Now suppose we are at 
a stage i € {1, ■ ■ ■ ,n— 2}. Again it holds that after embedding Vi+i the edges incident to Vi+i have to be 
consecutive. We call the set of leaves of Tj that correspond to these edges the black leaves. To go to the 
next stage, we thus first reduce T% so that all black edges appear together. A non-leaf node in the reduced 
PQ-tree is said to be black if all its descendants are black edges. We next create a new P-node Pi+\ and 
add all the out-edges of Vi + \ as its children. Now Tj+i is constructed from Tj as follows: 
Case 1: Tj contains a black node x that is an ancestor of all the black leaves. We obtain Tj+i from Tj by 
replacing x and all its descendants with p$+i. 

Case 2: T\ contains a (non-black) Q-node containing a (consecutive) sequence of black children. We obtain 
Tj+i from Tj by replacing these black children (and their descendants) with pj+i. 

This captures again exactly all possible circular orders of half-embedded edges around the newly em- 
bedded connected component. Note that as before the graph is non-planar if and only if the reduction step 
fails at any stage. 



4 Simultaneous planarity when the common graph is 2-connected 

This section contains our main result, a linear time algorithm for testing simultaneous planarity of two 
graphs that share a 2-connected subgraph. 

Let G\ = (V\,Ei) and G2 = (V 2 ,E 2 ) be two planar connected graphs with |Vi| = n\ and | V2 1 = n 2 . 
Note that G\\V\ n V 2 ] need not be the same as G 2 [Vi n V 2 ). Let G = (Vi D V 2 , E x n E 2 ) be 2-connected and 
n = \ V]_r\V 2 \. Let vi, v 2 , . . . , v n be an st-ordering of V\ Pi V 2 . We call the edges and vertices of G common 
and all other vertices and edges private. 

We say two linear or circular orderings of elements with some common elements are compatible if the 
common elements appear in the same relative order in both orderings. Similarly we say two combinatorial 
planar embeddings of Gi and G 2 respectively are compatible if for each common vertex the two circular 
orderings of edges incident to it are compatible. 

If G\ and G 2 have simultaneous planar embeddings, then clearly they have combinatorial planar em- 
beddings that are compatible with each other. The converse also turns out to be true, if the common 
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edges form a connected graph. This can be easily proved as follows. Let £\ and £2 be the compatible 
combinatorial planar embeddings of G\ and G2 respectively. Let £ p be the partial embedding of £\ [or £2] 
obtained by restricting G\ [resp. G2] to the common subgraph. (Note that since £\ and £2 are compatible, 
the partial embedding of £\ restricted to the common subgraph is the same as the partial embedding of 
£2 restricted to the common subgraph.) Now we can find a planar embedding of £ v and iteratively extend 
it to an embedding of £\ and an embedding of £2 (see Lemma 2 of Junger and Schultz [27] for a proof). 
The two planar embeddings thus obtained are simultaneous planar embeddings and thus it is enough to 
compute a pair of compatible combinatorial planar embeddings. 

We will find compatible combinatorial planar embeddings by adding vertices one by one. At any point 
we will have two sets of PQ-trees representing the partial planar embeddings of the subgraphs of G\ and 
G2 induced by the vertices added so far. Each PQ-tree represents one connected component of the current 
subgraph of G\ or G2. In the first phase we will add all private vertices of Gi and G2, and in the second 
phase we will add the common vertices in an st-order. When a common vertex is added, it will appear in 
two PQ-trees, one for G\ and one for G2 and we must take care to maintain compatibility. Note that if 
V\ = V2, i.e., when all vertices are common, only the second phase is needed and the whole algorithm will 
solely operate on two PQ-trees - one for each graph. 

Before describing the two phases, we give the main idea of maintaining compatibility between two 
PQ-trees. Recall that in Section |2.1| we found compatible orders for two PQ-trees using projection and 
intersection of PQ-trees. However, we were unable to store a set of compatible orderings as a PQ-tree, 
which is what we really need, since planarity testing involves building a sequence of PQ-trees as we add 
vertices one by one. 

To address this issue we introduce a boolean "orientation" variable attached to each Q-node to encode 
whether it is ordered forward or backward. Compatibility is captured by equations relating orientation 
variables. At the conclusion of the algorithm, it is a simple matter to see if the resulting set of Boolean 
equations has a solution. If it does, we use the solution to create compatible orderings of the Q-nodes of 
the two PQ-trees. Otherwise the graphs do not have simultaneous planar embeddings. 

In more detail, we create a Boolean orientation variable f(q) for each Q-node q, with the interpretation 
that f(q) = true if and only if q has a "forward" ordering. We record the initial ordering of each Q-node 
in order to distinguish "forward" from "backward". During PQ-tree operations, Q-nodes may merge, and 
during planarity testing, parts of PQ-trees may be deleted. We handle these modifications to Q-nodes 
by the simple expedient of having an orientation variable for each Q-node, and equating the variables as 
needed. When Q-nodes qi and q2 merge, we add the equation f(qi) = f{q2) if qi and q2 are merged in the 
same order (both forward or both backward), or f(qi) = ~~ , /((/2) otherwise. 

We now describe the two phases of our simultaneous planarity testing algorithm. To process the private 
vertices of G\ and G2 in the first phase we compute for each graph a reverse depth-first search ordering 
by contracting G into a single vertex and then running a depth-first search from this vertex. With these 
orderings we can now run the algorithm of Haeupler and Tarjan for all private vertices as described in 
Section [3j 

Now the processed vertices induce a collection of components, such that each component has an out- 
edge to a common vertex. Further, the planarity test provides us for each component with an associated 
PQ-tree representing all possible cyclic orderings of out-edges for that component. For each component 
we look at the out-edge that goes to the first common vertex in the st-order and re-root the PQ-tree for 
this component to have this edge represented by the special leaf. This completes the first phase. 

For the second phase we insert the common vertices in an st-order. The algorithm is similar to that 
described in Section[3jfor an st-order but, in addition, has to take care of merging in the private components 
as well. We first examine the procedure for a single graph. Adding the first common vertex v \ is a special 
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set-up phase which we describe first; we will describe the general addition below. 

Adding v% joins some of the private components into a new component G\ containing v\. For each 
of these private components we reduce the corresponding PQ-tree so that all the out-edges to v\ appear 
together, and then delete those edges. Note that due to the re-rooting at the end of the first phase 
the special leaf is among those edges. Thus the resulting PQ-tree represents the linear orderings of the 
remaining edges. We now build a PQ-tree representing the circular orderings around the new component 
C\ as follows: we take v\v n as the special leaf, create a new P-node as a root and add all the out-edges of 
v\ and the roots of the PQ-trees of the merged private components as children of the root (see Figure [I]). 

T ( C l) v \ v n 



T(C l ) TiC 1 ) T(C k ) 

Figure 1: Setting up T(Ci). The P-node's children are the outgoing edges of v\ and the PQ-trees for the 
components that are joined together by v\. 

Now consider the situation when we are about to add the common vertex Vi, i > 2. The graph so far 
may have many connected components but because of the choice of an st-ordering, all common vertices 
embedded so far are in one component Cj_i, which we call the main component. When we add Vi, all 
components with out-edges to V{ join together to form the new main component Cj. This includes Cj_i and 
possibly some private components. The other private components do not change, nor do their associated 
PQ-trees. 

We now describe how to update the PQ-tree Tj_i associated with C%—\ to form the PQ-tree T% associated 
with Cj. This is similar to the approach described in Section [3j We first reduce TJ-i so that all the black 
edges (the ones incident to v{) appear together. As before, we call a non-leaf node in the reduced PQ-tree 
black if all its descendants are black leaves. For any private component with an out-edge to Vi, we reduce 
the corresponding PQ-tree so that all the out-going edges to Vj appear together, and then delete those 
edges. We make all the roots of the resulting PQ-trees into children of a new P-node pi, and also add all 
the out-going edges of Vi as children of pi. It remains to add pi to Tj_i which we do as described below. 
In the process we also create a black tree Ji that represents the set of linear orderings of the black edges. 
Case 1: Tj_i contains a black node x such that all black edges are descendants of x. Let Ji be the subtree 
rooted at x. We obtain Ti from Tj_i by replacing x and all its descendants with pi. 

Case 2: TJ-i contains a non-black Q-node x that has a sequence of adjacent black children. We group 
all the black children of x and add them as children (in the same order) of a new Q-node x' . Let Ji be 
defined as the subtree rooted at x' . We add an equation relating the orientation variables of x and x' . We 
obtain Tj from Tj_i by replacing the sequence of black children of x (and their descendants) with pi (see 
Figure [2]). 

Note that we use orientation variables above for a purpose other than compatibility. (We are only 
working with one graph so far.) Standard planarity tests would simply keep track of the order of the 
deleted subtree Ji in relation to its parent. Since we have orientation variables anyway, we use them for 
this purpose. 
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Figure 2: (left) Adding vertex v% which is connected to main component Cj_i and to private components 
C 1 , . . . ,C k . (right) Creating Tj from Tj_i by replacing the black subtree by a P-node whose children are 
the outgoing edges of v i and the PQ-trees for the newly joined private components. 



We perform a similar procedure on graph G<i- We will distinguish the black trees of G\ and G2 using 
superscripts. Thus after adding Vi we have black trees J\ and Jf. It remains to deal with compatibility. 
We claim that it suffices to enforce compatibility between each pair J\ and Jf. 

To do so, we perform a unification step in which we add equations between orientation variables for 
Q-nodes in the two trees. 

Unification step for stage i 



We first project J\ and Jf to the common edges, as described in Section 2.1, carrying over orientation 
variables from each original node to its copy in the projection (if it exists). Next we create the PQ-tree Ri 
that is the intersection of these two projected trees as described in Section [2. 1| Initially Ri is equal to the 
first tree. The step dealing with Q-nodes (Step 3) is enhanced as follows: 

3. For each Q-node q of the second tree, and for each pair a\, a<i of adjacent children of q do the following: 
Reduce Ri by adding a consecutivity constraint on all the descendant leaves of a\ and 02- Find the 
Q-node that is the least common ancestor of the descendants of a\ and 02 in Ri. Add an equation 
relating the orientation variable of this ancestor with the orientation variable of q (using a negation 
if needed to match the orderings of the descendants) . 

Observe that any equations added during the unification step are necessary. Thus if the system of 
Boolean equations is inconsistent at the end of the algorithm, we conclude that G± and G2 do not have 
a compatible combinatorial planar embedding. Finally, if the system of Boolean equations has a solution, 
then we obtain compatible leaf-orders for each pair jj and Jf as follows: Pick an arbitrary solution to 
the system of Boolean equations. This fixes the truth values of all orientation variables and thus the 
orientations of all Q-nodes in all the trees. Subject to this, choose a leaf ordering / of Ri (by choosing the 
ordering of any P-nodes). I can then be lifted back to (compatible) leaf-orders of J\ and Jf that respect 
the ordering of /. The following lemma shows that this is sufficient to obtain compatible combinatorial 
planar embeddings of G\ and G2 

Lemma 1. If the system of Boolean equations has a solution then G\ andG2 have compatible combinatorial 
planar embeddings. 

Proof. The procedure described above produces compatible leaf orders for all pairs of black trees J/ and 
jf. Recall that the leaves of J\ (resp. Jf) are the out-edges of the component in G\ (resp. G2) 

and contain all the common in-edges of Vi. Focusing on G\ individually, its planarity test has succeeded, 
and we have a combinatorial planar embedding such that the ordering of edges around Vi contains the leaf 



9 



order of J*. Also, we have a combinatorial planar embedding of G2 such that the ordering of edges around 
Vi contains the leaf order of Jf. 

The embedding of a graph imposes an ordering of the out-edges around every main component. We 
can show inductively, starting from i = n, that the ordering of the out-edges around the main component 
Cj_i in G\ is compatible with the ordering of the out-edges in the corresponding main component in G<i. 
Moreover all the common edges incident to Vi belong to either Cj_i or Cj. This implies that in both 
embeddings, the orderings of edges around any common vertex are compatible. Therefore G\ and G2 have 
compatible combinatorial planar embeddings. □ 



4.1 Simultaneous planarity of k graphs 

In this subsection we consider a generalization of simultaneous planarity for k graphs, when each vertex 
[edge] is either present in all the graphs or present in exactly one of them. Let G\ = (Vi, G2 = 
(V2, E2), • • • , Gf. = (Vfc, Ek) be k planar graphs such that V = Vi n Vj for all distinct i,j and E = EiCi Ej 
for all distinct As before, we call the edges and vertices of G = (V, E) common and all other edges and 
vertices private. We show that the algorithm of Section [4] can be readily extended to solve this generalized 
version, when the common graph G is 2-connected. 

If G\, G2, ■ ■ ■ , Gk have simultaneous planar embeddings then they clearly have mutually compatible 
combinatorial planar embeddings. Conversely if Gi, G2, ■ ■ ■ , have combinatorial planar embeddings 
that are mutually compatible, then, as before (see the beginning of Section [4]), we can first find the planar 
embedding of the common subgraph and extend it to the planar embeddings of Gi, Gg, . . ■ , G^. Thus once 
again, the problem is equivalent to finding combinatorial planar embeddings for Gi, G2, . . . Gfc, that are 
mutually compatible. 

Our algorithm for finding such an embedding works as before, inserting private vertices first, followed by 
common vertices. The only difference comes in the unification step, where we have to take the intersection 
of k projected trees instead of 2. Doing this is straightforward: We initialize the intersection tree to be the 
first projected tree, and then insert the constraints of all the other trees into the intersection tree. Finally, 
Lemma [T] and its proof extend to multiple graphs. 



4.2 Running time 

We show that our algorithm can be implemented to run in linear time. (In the generalization to k graphs, 
the run time is linear in the sum of the sizes of the graphs — in other words, a common vertex counts k 
times.) Computing the reverse depth-first search ordering and the st-ordering are known to be doable 



in linear time 16 . The first phase of our algorithm uses PQ-tree based planarity testing with a reverse 
depth-first search order [22], which runs in linear time using the efficient PQ-tree implementation of Booth 
and Lueker [7j|8] . The re-rooting between the two phases needs to be done only once and can easily be done 
in linear time. The second phase of our algorithm uses PQ-tree based planarity testing with an st-order, 
as discussed in Section [3j This avoids re-rooting of PQ-trees, and thus also runs in linear time (§[22j[29 



The other part of the second phase is the unification step, which is only performed on the black trees, 
i.e. the edges connecting to the current vertex. Note that these edges will get deleted and will not appear 
in subsequent stages. Thus we can explicitly store the black trees and the intersection tree at every stage 
and allow the unification step to take time linear in the total size of both black trees. The intersection 



algorithm can be implemented in linear time, as mentioned in Section 2.1 The last thing that needs to be 
implemented efficiently is the handling of the orientation variables. It is easy to see that once the equations 
are generated, they can be solved in linear time, by repeatedly fixing the value of a free variable (true or 
false), finding all the equations that contain the variable and recursively (say in a depth-first way) fixing 
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all the variables so as to satisfy the equations. In the following sub-section we explain how to also generate 
the variable equations in linear time. With this we conclude that the algorithm runs in linear time. 

4.2.1 Generating variable equations in linear time 

Note that in the implementation of PQ-trees (see Booth and Lueker [§]) the children of a Q-node are stored 
in a doubly-linked list and only the leftmost and rightmost children have parent pointers. Thus when two 
Q-nodes, one a child of the other, merge, we may not know the variable and the orientation of the parent 
Q-node. To address this problem, we use labels on certain links of the doubly-linked list as explained below 
and compute all the equations generated by the reductions of a unification step at the end of the step. 

For any two adjacent child nodes Cj and a+i of a Q-node q, either the links Cj — > Cj+i and Cj+i — > Ci 
are labeled with / and —>l (respectively), for some literal I, or they are both unlabeled. The underlying 
interpretation is that q appears before Cj+i in the child ordering of q iff I is true. Thus the literals that 
we encounter when traveling from one end to the other of the doubly-linked list, are all (implicitly) equal. 
When a Q-node is first created with two child nodes, say x and y, we create a variable associated with it 
and label the link from x to y with the variable and the link from y to x with the negation of the variable. 

During the algorithm, there are two types of equations: (1) Equations consisting of literals appearing 
in Q-nodes of distinct trees. (These can be PQ-trees of the same graph, as happens in Case 2 of Section [4] 
or PQ-trees of different graphs, as happens in step 3 of Unification.) (2) Equations consisting of literals 
appearing in Q-nodes of a single tree (created during PQ-tree reductions). 

Note that type 1 equations are essentially equations that constrain the ordering of child nodes across 
the two Q-nodes, and we handle them as follows. Let c\,C2 be any two adjacent child nodes of the first 
Q-node that are constrained to appear in the same order as child nodes c' t , c' 2 of the second Q-node. If the 
links between c\ and C2 are unlabeled, we create a new variable x and label the links c\ — > C2 and C2 — > c\ 
with x and —>x respectively. Similarly, we label the links between d x and c 2 , if they are unlabeled. Now 
we create the equation that equates the literal associated with c\ — > C2 with the literal associated with 
c'i -»■ 4 

Type 2 equations happen when two Q-nodes merge. In this case the merged node contains literals 
from both the Q-nodes. Finding an equation after each merge is costly, as we need to scan through each 
Q-node until we find a literal. However we can compute the equations in a lazy fashion at the end of each 
unification step as follows. For every Q-node of the two black trees and the intersection tree obtained from 
their projections (i.e. the output tree of the unification step on the black trees), we pass from the first child 
to the last child and equate all the literals encountered in the labels of the links. This clearly takes linear 
time in the size of the black trees. 

4.3 Simultaneous planar drawings 

Our algorithm to test simultaneous planarity finds compatible combinatorial planar embeddings but does 
not find actual simultaneous drawings. In this section we make some remarks on finding simultaneous 
planar drawings. 

Two natural goals are to minimize the number of bends along the polygonal curves representing the 
edges, and to minimize the number of crossings between the private edges of the two graphs. We will show 
that we can find drawings where each edge has at most n bends, and each pair of private edges intersects 
at most once. Recall that n = \V± PI V2I. In fact any edge of G\, and any private edge of G2 that joins 
private vertices of G2 will be drawn as straight line segments. 

Take a straight-line planar drawing of the first graph, G\. This includes a straight-line drawing of the 
common graph G, which we will extend to a planar drawing of G*2. 
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The main idea is to draw the private parts of G2 inside the faces of the drawing of G using shortest 
paths. (Of course, we must offset the shortest paths slightly to avoid coincident features.) A shortest path 
inside a polygon has two crucial properties: it bends only at reflex vertices of the polygon, and it intersects 
any line segment contained in the polygon at most once. 

Any bend on an edge of G2 will occur in the neighborhood of a common vertex (a vertex of G). We 
will maintain the following invariants: (1) each edge of G2 has at most one bend per common vertex; (2) 
the boundary of each face of G2 has at most one reflex bend per common vertex; and (3) any edge of G\ 
intersects any face of G2 in a single line segment. Property (1) implies that each edge of G2 has at most n 
bends. Property (2) allows us to maintain property (1) as we draw more of G2. Property (3) ensures that 
private edges of Gi cross private edges of G2 at most once. Initially, the only part of G2 that is drawn is 
G, and the properties are satisfied. 




Figure 3: (left) A face of G (drawn with solid lines), some private edges of G\ (dotted lines), and two 
private edges of G2 (dashed lines) with endpoints in G: (d,a) and (c, g), each drawn with a bend in the 
neighborhood of the common vertex e. Observe that although there are two bends near e, each face of G2 
has only one reflex bend near e. (right) Adding private vertices of G2: vertex v is drawn along a shortest 
path from a to b and a component is drawn in a small disc around v. 

We first draw the private edges of G2 whose endpoints are in G. Consider one such edge. The 
combinatorial embedding specifies the face, /, of the partial embedding of G2 that contains the edge. Face 
/ is drawn as a polygon and, by property (3), any reflex vertex of the polygon is in the neighborhood of a 
unique common vertex. Draw the edge as a shortest path in /, pulling it slightly away from the boundary 
of /. Observe that this satisfies our invariants. 

We now insert private vertices of G2 and their incident edges. Consider a connected component C of 
G2 — (V\ n V2), consisting of some private vertices and their induced edges (call this C) together with some 
edges joining C to the common graph G. The combinatorial embedding specifies the face, /, of the partial 
embedding of G2 that contains G. Face / is drawn as a polygon and, by property (3), any reflex vertex of 
the polygon is in the neighborhood of a unique common vertex. We will contract C to a single vertex v 
and draw v (as a point) together with all the edges joining it to common vertices. After that, we take a 
small disc around the point, small enough to avoid all edges of G\, and we make a straight line drawing of 
C inside the disc. Thus it suffices to consider adding a single private vertex v and the edges joining it to 
common vertices. If v is adjacent to only one common vertex, we can augment G2 to make it adjacent to 
another common vertex. Let a and b be two common vertices adjacent to v. Draw a shortest path inside 
/ from a to b, and choose some point along that path to place v. Draw the remaining edges incident to 
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v as shortest paths in the polygon, pulling the paths slightly away from the boundary of / and from each 
other to avoid coincident features. Observe that v is not a reflex vertex in any of the resulting faces. This 
construction satisfies our invariants: each edge has at most one bend per reflex vertex of /, the boundary 
of each resulting face has at most one reflex bend per vertex of /, and any edge of G\ intersects any face 
of Gi in a single line segment. 



5 Conclusions 

We have given a linear-time algorithm to test simultaneous planarity of two graphs whose common subgraph 
is 2-connected. Our algorithm does not require the two graphs to have the same vertex set. Furthermore, 
our algorithm works for the more general case when there are k graphs, any two of which have the same 
2-connected subgraph in common. 

We conjecture that simultaneous planarity of multiple graphs can be tested in polynomial time when 
any vertex/edge is either private to a single graph or common to all graphs. Our algorithm solves the 
special case when the common graph is 2-connected. Note that simultaneous planarity is known to be 
NP-hard for three graphs in general |20|. The three graphs constructed in the reduction have edges in all 
possible combinations: edges private to each graph, edges common to all three graphs, and edges common 
to a pair of graphs but not the third, for every pair. 

Another interesting case of simultaneous planarity for multiple graphs is when the graphs intersect in 
layers: the graphs are ordered G%, . . . , Gk, and every vertex/edge occurs in a consecutive subsequence of 
the graphs. This layered structure arises from a graph changing over time, so long as a vertex/edge does 
not re-appear after it vanishes. The complexity of layered simultaneous planarity is open, even for three 
graphs. For three graphs the layered structure is equivalent to the the condition that if a vertex/edge is in 
G± and G3 then it is also in G<i- 

A weaker version of simultaneous planarity for two graphs requires only that each common vertex be 
represented by the same point — a common edge may be represented by different curves in the two planar 
drawings. Any pair of planar graphs can be represented this way: after choosing arbitrary vertex positions, 
each graph can be drawn independently, as shown by Pach and Wenger |33| . It would be interesting to 
generalize our algorithm to the case where some designated common edges are allowed to be represented 
by different curves in the two drawings. There is also a natural optimization version of the problem: given 
two planar graphs with some common vertices and edges, find planar drawings of the graphs so that every 
common vertex is drawn as the same point in the two drawings, and maximize the number of common 
edges that are drawn as the same curve in the two drawings. 

Lastly, both SPQR-trees [5] and PQ-trees have been used for many planarity related problems and both 
give rise to distinct representations of all planar embeddings of a (2-connected) planar graph [5 11 22 



Understanding the strength of both representations and how they relate better is an important question. 
Comparing the PQ-tree approach to simultaneous planarity taken here with the SPQR-tree approach in (3] 
or [17] might be a good start. The recent manuscript |6| makes very interesting progress in this direction. 
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